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ABSTRACT 

This  paper  gives  a  mathematical  theory  for  a  visco-elastic  fluid 
exhibiting  a  hysterisis  loop  in  the  shear  stress  vs.  shear  rate  plane.  The 
main  rheological  idea  is  to  introduce  a  constitutive  equation  of  rate  type 
whose  steady  shear  stress  vs.  shear  rate  locus  is  non-monotone.  The  main 
mathematical  idea  is  to  use  a  local  shock  structure  theory  to  pick  out  the 
admissible  solutions  in  loading  and  unloading  of  the  applied  driving  force. 
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SIGNIFICANCE  AND  EXPLANATION 


In  processing  polymers  a  main  concern  is  loss  of  stability  in  what  is 
often  termed  melt  fracture.  This  loss  of  stability  is  often  thought  of  as 
being  due  to  internal  material  properties.  This  paper  analyzes  how  a  certain 
material  model  for  a  visco-elastic  fluid  yields  results  with  some  features  of 
the  "ripple"  phenomena  seen  in  melt  fracture  of  certain  polymers. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 


UNSTABLE  VISCO-ELASTIC  FLUID  FLOW 
EXHIBITING  HY&tERITIC  PHASE  CHANGES 

•  ** 

J.  K.  Hunter  and  M.  Slemrod 

0.  Introduction 

In  his  1969  survey  paper  [1]  on  melt  fracture  of  molten  polymers 

J.  P.  Tordella  described  a  class  of  experiments  illustrating  a  phenomenon  he 

called  "ripple".  In  these  experiments  certain  polymers  (high  density 

polyethylene  and  crystalline  copolymers  of  tetrafluorethylene  and 

hexofluorophylene)  exhibit  remarkable  behavior  in  shear  flow.  While  this 

topic  has  been  discussed  many  places  in  the  rheological  literature  (see  [2] 

for  a  bibliography)  it  is  only  recently  [3] ,  [4]  that  mathematical  analysis 

has  been  used  in  the  qualitative  study  of  the  "ripple"  phenomenon.  In  this 

paper  we  continue  this  program  by  giving  a  solution  of  a  boundary  value 

problem  which  has  properties  similar  to  the  behavior  described  in  [ 1 ] . 

First  we  paraphrase  Tordella' s  observations.  For  shear  flow  in  a 

capillary  tube  at  sufficiently  small  values  of  stress  ordinary  visco-elastic 

flow  occurs:  viscosity  decreases  with  increasing  stress  and  the  emerging 

stream  is  smooth  and  regular.  At  and  above  a  critical  stress  T  the 

A 

emerging  stream  has  some  distortion.  At  a  second  critical  stress  T  which 

B 

is  greater  than  x  the  shear  rate  tends  to  be  double  valued.  Furthermore 
A 
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upon  increasing  shear  stress  the  shear  rate  possesses  values  on  different 
branch  of  the  stress  vs.  rate  of  shear  curve*  Decreasing  stress  brings  a 
return  to  the  original  branch  in  a  hysteretic  fashion. 

In  this  paper  we  attempt  to  give  a  qualitative  explanation  for  this 
hysteresis  effect  based  on  an  elementary  visco-elastic  model  of  rate  type. 

The  main  feature  of  the  model  is  that  for  steady  flows  it  possesses  a  non- 
monotone  shear  stress  vs.  rate  of  shear  constitutive  relation.  Such  models 
have  been  suggested  for  phase  transition-like  behavior  in  non-linear 
elasticity  by  Ericksen  [5]  and  more  recently  in  non-linear  visco-elasticity  by 
Bernstein  and  Zapas  [6] .  Within  the  range  of  visco-elastic  fluids  such  ideas 
seem  mainly  found  in  the  Soviet  literature,  i.e.  in  the  work  of  Vinogradov, 
Rutkevich  and  their  co-workers  [7],  [8],  [9],  [10].  What  makes  our  work 
different  is  we  solve  a  specific  boundary  value  problem  for  the  steady  shear 
flow  of  a  visco-elastic  fluid  satisfying  our  constitutive  relation.  The 
boundary  value  problem  (in  general)  possesses  non-unique  solutions  and  we  use 
a  local  dynamic  shock  structure  theory  to  choose  a  solution.  This  solution 
exhibits  a  hysteresis  loop,  double  valued  shear  rates  at  critical  stresses, 
and  spatially  segregated  flow  regimes  similar  to  the  flow  birefrigence 


phot graphs  given  in  [1] 


1 •  General  relations 

If  In  a  fixed  Cartesian  coordinate  system  x,  y,  z  the  velocity  field  of 
a  flowing  and  incompressible  fluid  body  has  the  form 

vx  -  0,  «  v(x,t),  v*  ■  0  ,  (1.1) 

we  say  the  motion  is  a  rectilinear  shearing  flow.  For  such  a  flow  the 
condition  of  incompressibility  div  _v  =  0  is  automatically  satisfied,  if  the 
fluid  is  a  simple  fluid  in  the  sense  of  Coleman  and  Noll  [11]  then  the 
components  of  the  Cauchy  stress  obey  the  relations 

T  =  -pi  +  T 

■v  *** 

where  the  extra  stress  T  satisfies 

*o  <V8))  ' 
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Here  p  is  an  indeterminate  hydrostatic  pressure  and  A q ,  4 1 ,  A  2  are  real 
valued,  generally  non-linear  functions  of  the  relative  shearing  history 


At(s), 

Afc(s)  =  -  /£_g  vx(x,T  )dr  (0  <  s  <  »)  . 

From  the  fluid's  isotropy  we  know  the  functionals  Aq,  a.,  A 2  satisfy  the 


relations 


(1.3) 


<-At(s))  -  -4q  ( Afc ( 8 ) )  , 

S*0  8*0 

"i  (-At(s))  *  "i  ( At < s ) ) ,  i  *  1,2  . 

8*0  8*0 


If  we  substitute  conditions  (1.2)  into  the  equations  of  balance  of  linear 
momentum  we  find 
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xx  yx 

where  b^  is  a  component  of  the  body  force  b  *  (bj,0,0)  and  Tg  ,  are 

functions  of  x  and  t  only.  From  (1.4a)  we  find  *  Y(t)  and  b.|(x,t) 
is  the  body  force  component  necessary  to  preserve  (1.4a)  as  an  identity. 

Thus  v(x,t)  satisfies  the  evolution  equation  (1.4b). 

To  proceed  further  it  is  necessary  to  make  some  mathematical  assumptions 
as  to  nature  of  the  functional  6q.  For  this  analysis  we  assume  <SQ  has  the 
particularly  simple  form 

T^(t)  -  “0  (At(s))  -  f~  e'Xs0(vx(X|t-s))ds  .  (1.5) 

s*0 


Alternatively  we  note  that  it  is  equivalent  to  (1.5)  that  T*y(t)  be  of  rate 

E 

xy 

type  since  T  satisfies 

E 

§£  T^(t)  +  XT^(t)  -  o(vx(*,t))  .  (1.6) 

In  the  case  y  is  a  constant  (representing  a  constant  applied  during 
force)  (1.6)  and  (1.4b)  imply  v(x,t)  satisfies  the  second  order  equation 


where  Y  “  Xy«  If  we  set  vfc  =  u,  vx  -  w  we  see  that  u,  w  satisfy  the 


first  order  system 


p(u  +  Xu)  -  o(w)  -  Y  i 


w^  =  u 
t  X 


We  will  assume  our  visco-elastic  fluid  is  confined  between  two  parallel 
walls  of  infinite  extent  at  x  =  -1  and  x  -  1  and  satisfies  no-slip 


boundary  conditions 


v  -  0  at  x  *  -1#  1  . 


(Of  course  this  implies  u  •  0  at  x  *  -1,  1.) 

So  far  we  have  made  no  assumptions  as  to  the  nature  of  a  beyond  the 
fact  that  (1.3)  forces  o  to  be  an  odd  function  of  its  arguments.  In  this 
paper  we  take  o  to  be  a  smooth  odd  function  with  graph  shown  in  Figure  1. 


i  a-phase 


4—  6-phase  ■ 


Here  o  is  such  that 

o'  >  0  on  [0,wa),  (Wg»“)  i 
o'  <  0  on  »wg)  , 

(1.10) 

°<V  -  '"'•s’  * T.  - 

0<w  )  -  o(w  )  -  Yk  . 
a  2d 

Following  the  standard  nomenclature  of  phase  transitions  we  say  that  at  a 

point  where  w  *  vx  takes  on  values  in  [  0 ,  w^  )  the  fluid  is  in  the  a -phase 

and  when  w  takes  on  values  in  (wa )  the  fluid  is  in  the  3-phase.  Of 

p 

course  for  different  values  of  x  the  same  fluid  may  be  in  the  a  and  3 
phases  simultaneously. 

Since  the  function  a  is  not  globally  invertible  we  denote  by  o”\ 

Oa  1  the  respective  inverses  of  a  in  the  a  and  3  phases  i.e. 

P 

a’V)  €  £°*wal  for  0  <  y  <  Yb  , 

0“\y)  €  [wp#“)  for  Ya  <  Y  <  -  . 

In  what  follows  we  admit  solutions  with  values  in  a  and  3  phases. 

However  we  shall  not  allow  solutions  with  values  in  (w  ,w. )  as  (1.7)  is  an 

a  3 

elliptic  initial  value  problem  in  this  region  and  hence  will  exhibit  the 
classical  Hadamard  instability. 


2*  Bifurcation  and  coalescence  of  fluid  phases  for  steady  motion 


In  this  section  we  will  show  how  discontinuities  between  a  and  0 
phases  may  appear  and  disappear  in  nearly  steady  flow  ( vfc  =  0).  For  steady 
flow  (1.7)  implies  0(v  )  «  Y  which  upon  integration  yields  <J(v  )  =  Yx  + 

X  X 

constant.  By  symmetry  of  the  flow  x  about  x  =  0,  v(x)  should  be  an  even 

function  of  x  and  hence  vx  should  be  odd.  As  noted  earlier  o  is  an  odd 

function  by  isotropy.  Hence  otv^)  -  Y*  is  odd  in  x.  The  only  constant 

which  is  odd  in  x  is  zero  so  at  equilibrium  (1.7)  reduces  to 

<j(v  )  «  Yx  .  (2.1) 

x 

If  0  <  Y  <  Y  then  0  <  Yx  <  Y  for  x  e  [0,1].  In  this  case  (2.1) 

S  cl 

has  a  unique  solution 


v(x)  =  J*  0a1(Ys)ds  , 
vx(t)  *  aa1(Yx)  . 

On  the  other  hand  for  Y  <  Y,  the  lack  of  a  unique  inverse  allows 

a 

to  possess  non-unique  solutions*  For  example  we  may  have 

* 

v(x)  -  /X#  o^1 (Ys)ds  +  /*  Og1(Ys)ds  , 
x 

-1  * 
vx«x)  -  0  (Yx)  ,  0  <  x  <  x  ; 

v ( x)  ■  /*  0^1(Y8)ds  , 


(2.2) 


(2.1) 


(2.3) 


v  (x)  -  ©I1 (Yx)  , 
x  p 


* 

X  <  x  <  1  . 


The  function  v(x)  given  by  (2.3)  is  continuously  differentiable  on  (0,x*), 
( x* , 1 ) ,  continuous  on  (0,1),  and  satisfies  the  boundary  conditions 
v(1)  m  0.  Of  course  the  values  for  -1  <  x  <  0  are  obtained  by  symmetry. 

As  x*  is  an  arbitrary  number  between  Y  /Y  and  1,  (2.3)  represents  an 

a 

infinite  number  of  solutions  of  (2.1)  possessing  a  jump  discontinuity  in  v 


at  x  .  We  note  that  on  (0,x  )  the  fluid  is  in  the  a -phase  and  on  (x  , 

* 

the  fluid  is  in  the  $ -phase.  So  the  line  x  =  x  represent  a  singular 
surface  where  the  fluid  exhibits  an  inter-phase  jump.  In  addition  to 
solutions  of  the  form  (2.3)  if  Y  <  Y  <  Y.  ,  (2.2)  will  provide  a 

d  D 

continuously  differentiable  solution  of  (2.1).  Pictorially  these  solutions 
are  illustrated  in  Figure  2. 


1) 


Figure  2 

(A)  Solution  (2.2),  0  <  Y  <  Y,; 

b 

(B)  Solution  (2.3),  Y  <  Y  . 

£1 


Finally  for  Y  >  Y^  (2.1)  can  only  possess  solution  existing  in  both  a 
and  $  phases.  In  this  a  jump  discontinuity  in  vx  at  some  inter-phase 
surface  is  inevitable. 

Now  consider  the  following  "quasi-static"  ejqperiment.  Let  Y  be  slowly 

raised  from  zero.  For  Y  €  [0,Y  ),  (2.1)  possesses  only  C1  solution  given 

£1 

by  ( 2. 2) .  When  Y  is  raised  above  Y  we  face  the  possibility  of  a  spatial 

d 

bifurcation  occuring  with  the  appearance  of  a  solution  of  the  form  (2.3). 
Finally  continuing  to  raise  Y  beyond  Y^  to  Yq  guarantees  the  formation 
of  a  solution  possessing  a  jump  discontinuity  in  vx.  If  we  now  decrease  Y 
from  Y  Q  back  to  zero  we  sh-  J.d  eventually  see  the  a -phase  solution  (2.2) 
when  Y  is  less  than  Yfl«  Thus  we  expect  to  see  bifurcation  and  coalescence 


3*  Local  dynamic  analysis 


Our  point  of  view  is  simple.  We  believe  that  for  a  spatial 
bifurcation/coalescence  of  fluid  phases  to  occur  near  equilibrium  an  inter¬ 
phase  singular  surface  across  which  vx  jumps  must  slowly  propagate  in  the 
fluid.  Hence  if  we  understand  when  such  surfaces  can  propagate  we  will  know 
when  a  spatial  bifurcation/coalescence  of  fluid  phases  can  take  place. 

We  shall  call  a  C1  surface  S  :  x  =*  s(t)  a  singular  surface  for  (1.8) 
if  u  and  w  experience  jumps  across  S.  If  S  is  a  singular  surface  let 
(8(t)/t)  be  a  fixed  point  on  the  graph  of  S  in  the  x  -  t  plane  and 
T  ■  s(t).  Denote  by  u+,  w+,  u_,  w_  the  respective  limits  fr-m  the  right  and 
left  as  (x,t)  -*■  (s(t),t)  for  u,  w  in  the  x  -  t  plane.  If  we  put 
[u]  “  u+  -  u_,  [w]  =  w+  -  w_,  etc.  then  we  know  the  Rankine-Huginot 
conditions  must  be  satisfied,  i.e. 


pT[u]  +  [a]  -  0 
T[w]  +  [u]  -  0 


(3.1) 


Of  course  (3.1)  implies 

pt2“M*  (3*2) 

Our  goal  is  to  find  out  which  singular  surfaces  for  (1.8)  are  physically 

meaningful.  To  do  this  we  adopt  the  following  premise:  The  shearing  stress 

given  by  (1.5)  possesses  no  Newtonian  viscous  contribution.  We  believe, 

however,  that  any  real  fluid  should  possess  some  Newtonian  contribution  albeit 

small.  Hence  the  "good"  solutions  (1.7)  and  (1.8)  should  be  the  ones  that  are 

robust  or  stable  (in  sense  to  be  made  precise  below)  under  small  perturbations 

of  (1.7)  and  (1.8)  with  a  Newtonian  viscous  term.  Of  course  this  view  is  not 


original  and  has  a  long  tradition  dating  back  to  Rayleigh  [12].  More  recent 


discussions  of  the  role  of  viscosity  in  the  study  of  singular  surfaces  may  be 
found  in  the  books  of  Zel'dovich  and  Raizer  [13]  and  Courant  and  Friedrichs 
[14]. 

If  we  add  a  linear  Newtonian  viscous  term  yvx  to  the  shearing  stress 
with  y  >  0  denoting  an  (assumed)  small  viscosity  we  find  (1.5)  is  replaced 
by 

T^tt)  =  /n  e  *so(v  (x,t-s))ds  +  yv  .  (3.2) 

Substituting  (3.2)  into  the  balance  of  linear  momentum  equation  yields  the 
evolution  equation 

pvt  =  r0  e"Xso(vp<x,t-s)),xds  +  yvPx  -  y  (3.3) 

where  vP  is  used  to  denote  the  solution  of  (3.3).  If  we  set  wP  =  vP, 

u‘  =  vP,  we  readily  see  by  differenting  (3.3)  with  respect  to  t  and 

y  y 

integrating  by  parts  that  w  ,  u  satisfy  the  first  order  system 
p(uH  +  Au**)  =  <XU(wU)  +  yuW  -  Y  »  -1  <  x  <  1 

t  X  XX 

(3.4) 


where  op(w)  is  defined  by 

<JP(w)  «=  0(w)  +  Ayw  . 

We  now  follow  an  argument  of  Dafermos  [15]  for  the  analysis  of  singular 
surfaces.  Let  S  be  a ‘singular  surface  for  (1.8)  and  (s(t),t),  u+,  w+,  u_, 
w_,  T  as  above.  As  we  are  interested  in  the  behavior  of  u,  w  in  the 
immediate  behavior  of  S  we  introduce  stretched  normal  and  tangential 
coordinates 

5  -  x-s(t)-T( t-t)  n  *  (T(x-s(t))+(t-t) 

/  _  • 

V,  M 

Let  u^tC^n)  *  uW(x,t),  ^(£,0)  ■  wP(x,t).  Then  an  application  of  the  chain 


rule  shows  that  (3.4)  is  equivalent  to 


1 


I 


3 


p(-Tu*i  +  u11  +  Xyu^)  * 

t  n 

o“V)(^  +  w^T)  -  Xy(^  +  TWj) 

(3.5) 

t  (3J5  *  zr^n  ♦  T^n)  -  „r 

-tv*  ♦  ?  -  S  .  tS”  . 

€  n  £  n 


A  classical  result  of  Maxwell  [16]  asserts  that  the  discontinuities  in  u 

and  w  can  occur  only  normal  to  the  singular  surface  S.  Hence  if  as 

y  ♦  0+  we  desire  u^,  wW  to  approach  the  discontinuous  profile  of  u,  w, 

-y 

u  ,  w  should  tend  to  functions  which  change  only  in  a  direction  normal  to 

the  singular  surface.  In  other  words  we  expect  uw,  ww  ♦  u®,  w®  as  y  +  0+ 
— 0  — 0  . 

where  u  ,  w  are  functions  of  £  alone*  Thus  formally  taking  limits  as 
y  ♦  0+  in  (3.5)  we  find 


-Tp  u°(5)  *  0(w°)*  +  u°(C) 
-T  ^°(«)  -  ?(«  , 


(3.6] 


d  y  y 

where  •  ■  in  addition  for  u  ,  w  to  approximate  the  discontinuous 

profile  of  u,  w  we  will  need  lim  uFtSrO)  -  u  ,  lim  u^^O)  *  u  , 

5+0-  "  5+0+ 

lim  w^tCfO)  *  w  ,  lim  ( C # 0 )  m  w  .  If  we  now  formally  take  the  limit  as 

5+0-  ”  5+0+ 

y  ♦  0+,  interchanging  limiting  processes,  we  find  u  ,  w  should  satisfy 
boundary  conditions 


— 0 

u  (-  •)  -  u  ,  w  (-  •)  «  w_  , 

— 0  wj 

W  (+  •)  ■  u+,  w  (+  *)  -  w+  . 

This  discussion  motivates  the  following  definition. 


(3.7) 


2- 


Definition.  A  singular  surface  S  :  x  *  s(t)  for  (1.8)  with  (s(t),t),  u_, 
w_,  u+,  w+,  T  as  above  is  admissible  according  to  the  Newtonian  viscosity 
criterion  if  (3.6)  possesses  a  C2  solution  for  every  point  (s(t),t)  on  the 
graph  S. 


It  will  be  the  ability  to  satisfy  the  Newtonian  viscosity  criterion  that 
will  characterize  the  robustness  or  stability  of  the  shearing  stress  under 
small  perturbations  of  Newtonian  viscous  stress. 

Theorem  3.1.  As  singular  surface  S  for  (1.8)  is  admissible  according  to  the 
Newtonian  viscosity  criterion  if  and  only  if 


o(w)  -  o(w  ) 

- s - —  -  PT2 

w-w 


>  0 

if 

T  > 

0 

<  0 

if 

T  < 

0 

(3.8) 


for  every  w  between  w_  and  w+.  Alternatively  if  (w+-w_)T  >  0  (or 
(w+-w_)T  <  0)  the  chord  which  joins  (w_,o(w_))  to  (w+,o(w+))  lies  below 
(above)  the  graph  of  a  between  w_  and  w+. 

Proof .  The  proof  has  been  given  in  [15]  and  [17].  However  since  it  is  short 
and  straightforward  we  present  it  for  completeness.  First  integrate  (3.6) 
from  -  00  to  C •  We  see  that  w^(C)  satisfies  the  first  order  ordinary 
differential  equation 

T2p(w°(£)  -  w_)  -  0(w°)  +  0(w_)  +  TW° ( K)  ■  0  (3.9) 

•“0 

with  boundary  conditions  w  (-  •)  *  w_,  w  (+  •)  ■  w+.  If  w_  <  w+  a 

necessary  and  sufficient  for  there  to  exist  an  orbit  connecting  w_  to  w+ 

2  *  *  * 
is  that  T{T  p(w  -  w_)  -  o(w)  +  o_)  <  0  for  w_  <  w  <  w+  for  this  gives  a 

one  directional  vector  field  between  w_  and  w+.  If  w_  >  w+  the  sign 

changes  in  the  inequality.  But  this  is  just  a  restatement  of  (3.8). 

Corollary  3.2.  If  w_  is  in  the  a-phase  and  w+  is  in  the  0-phase 


with  w_  <  w+  then  (3.8a)  ((3.8b))  holding  means  the  singular  surface 


propagates  into  the  3 -phase  (a -phase)  with  the  a-phase  (3 -phase)  behind 
the  singular  surface. 

Proof.  This  follows  from  the  definition  of  C * 

Mow  let  us  reconsider  the  "quasi-static"  experiment  in  light  of  our 
admissibility  criterion.  We  have  seen  that  as  we  slowly,  monotonically  raise 
Y  from  zero  beyond  y  ,  but  less  than  y.  ,  there  exists  the  possibility  of 

A  D 

(2.1)  possessing  both  a  single  a-phase  solution  (2.2)  or  two  phase  solutions 
e.g.  of  the  type  (2.3).  But  for  a  two  phase  solution  to  form  an  inter-phase 
singular  surface  between  a-  and  3 -phases  must  appear  in  the  original 
strictly  a-phase  solution  with  the  3 -phase  propagating  into  the  a-phase. 

If  y  is  raised  slowly  this  inter-phase  singular  surface  is  expected  to 
propagate  slowly.  So  we  are  led  to  ask:  when  can  we  have  an  admissible 
slowly  moving  singular  surface  separating  a-  and  3 -phases  with  3 -phase 
propagating  into  a-phase.  According  to  Corollary  3.2  a  slowly  moving  surface 
of  this  type  must  possess  left  and  right  limits  (in  the  x-t  plane)  near  wq 
and  W2>  See  Figure  3 


Figure  3 

3 -phase  slowly  into  a-phase 
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What  is  the  location  of  the  inter-phase  singular  surface?  If  the  surface 

were  actually  static  with  T  =  0  and  located  at  x  *  x*  >  0  with  w 

possessing  limits  wq  and  w2  as  xS*x*  and  x^x*  then  (2.1)  implies 
-1  *  -1  * 

w  *  o  (yx  )  and  w.  *  ae  (yx  ).  Since  o(w  )  =  c( w„)  =*  Y.  we  find 
a  a  2  p  a  2  b 

* 

x  *  Y./Y*  Using  the  sytmnetry  of  v  with  respect  to  x  we  then  see  that 
x> 

admissible  inter-phase  singular  surfaces  form  when  y  =*  Y^  at  x  =  ±1  and 
propagate  with  $-phase  moving  into  the  channel.  This  is  shown  in  Figure 
4.  Of  course  if  Y  increases  the  location  of  x*  asymptotically 
approaches  x  -  0.  Also  since  T  is  not  zero  but  merely  small  and  negative 
x*  represents  only  an  approximation  to  the  location  of  the  singular  surface. 


/ 

s 

x  =  1 

B-phase 

x  =  x*  =  Yb/Y 

a-phase 

J,T<  0 

x  =  0 

X  =  -X*  =  -Yb/Y  \ 

g-phase 

x  =  -1 

Figure  4 

Solution  of  (2.1)  for  Y  increasing,  Y  >  Y.  • 

D 
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We  now  move  onto  a  second  question:  when  can  we  have  a  slowly  moving 


inter-phase  singular  surface  with  a-phase  moving  into  3-phase?  By  an 
argument  analogous  to  the  one  given  above  we  see  this  can  only  occur  when  the 
limits  from  the  left  and  right  are  near  w^  and  w^.  This  is  illustrated  in 
Figure  5. 


Figure  5 


What  is  the  location  of  this  inter-phase  singular  surface?  By  an 

argument  analogous  to  the  one  given  above  we  find  it  is  approximately  located 
** 

at  x  -  Y_/Y«  We  can  see  the  effect  of  slowly  monotonically  lowering  Y 

from  some  value  Yq»  Yq  >  Y^»  in  our  "quasi-static"  experiment.  The  singular 

* 


surface  is  initally  located  at  x  =»  YK/Yn  (say  as  shown  in  Figure  4  with 


Y  “  Y0)»  As  we  lower  Y  the  inter-phase  singular  surface  cannot  move  until 


Vo  Vo 

— —  The  reason  is  that  when  — —  <  Y  <  Yn  if  the  surface  (for  x  >  0) 


'b  'b 

is  to  move  upward  it  must  occur  with  1  >>  T  >  0  and  a-phase  moving  into 


3-phase.  But  we  have  just  observed  this  can  only  occur  at  x  *  Y  /Y« 

A 


**  *  ^Vo 

However  x  *  xQ  can  only  happen  when  y  =»  — - —  is  a  contradiction.  On  the 

b 
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other  hand  if  the  surface  (for  x  >  0)  is  to  move  downward  we  must  have 

6-phase  propagating  into  a-phase  with  T  <  0  and  small.  This  of  course 

* 

must  take  place  at  x  *  Y^A  with  increasing  y  not  decreasing  y.  So  this 

* 

possibility  is  also  excluded.  So  the  surface  remains  at  x^  -  Y^/Yq  until 

Vo 

Y  reaches  - •  Thus  the  solution  exhibits  a  form  of  shape  memory  in  that 

'b 

interface  "remembers"  the  location  it  was  brought  to  via  the  Y  increasing 
part  of  the  cycle. 


As  we  bring  Y  below 


the  singular  surface  now  begins  to  move 


slowly  upward  with  a-phase  propagating  into  6-phase.  The  surface  is 

** 

approximately  located  at  x  =  Y  A«  When  Y  reaches  Y  the  inter-phase 

A  ft 

singular  surface  moves  into  the  channel  walls  and  disappears.  So  for 
0  <  Y  <  Y  we  see  only  the  a-phase  solution  ( 2. 2) .  Our  admissible  solution 

ft 

to  (2.1)  is  summarized  in  the  following  table. 


Table: 

Solution  of  (2.1)  for  Y  increasing  monotonically  from 


0  to  yq  >  Y^#  then  decreasing  monotonically  to  zero 
( 1 )  Y  increasing,  0  <  Y  <  Y.  : 

D 


V(x)  -  /*  Oa1(Ys)ds,  0  <  X  <  1  , 


v\(x)  =  a_  (yx)  i 


(2)  Y  increasing,  y  <  y  <  YQ  s 


Y  /Y 

v(x)  -  J*  Oa1(Ys)ds  +  Cg'tYsJds, 

b 

v  (x)  -  (T^yx)  , 
x  a 

y 

v(x)  =  J*  Og1  (Ys)ds,  <  x  <  1  , 

v  (x)  =»  a o^Yx)  ; 
x  p 

(3)  v  decreasing,  —  Y0  <  Y  <  Yn  : 

Tb 

Y  /Y 

v(x)  =  /*  .  o  ^  (Ys)ds  +  /  b  °  O-^Yslds, 

WT0  p 

vx(x)  -  o'^Yx)  , 
v(x)  -  f*  Og1 (Ys)ds, 
v  (x)  *  o'^ox)  ; 

X  p 

(4)  y  decreasing,  Yfl  <  Y  <  ~  YQ  « 

b 

Y  /Y 

V(x)  =  /*  /Y  o^1  (Ys) ds  +  Og^YsJds, 

d 

vx(x)  “  °a1(Yx>  ' 

v{x)  »  /*  Og 1 (Ys)ds,  Ya/Y  <  x  4  1  , 


><.«i  . 
T0 


4.  Resolution  of  a  inter-phase  discontinuity  in  the  neighborhood  of  the 

singular  surface 

In  the  previous  section  we  claimed  that  for  w  values  in  the  a-phase 
near  wq  and  w  values  in  the  3-phase  near  Wj  there  will  be  a  slowly 
moving  inter-phase  singular  surface  propagating  from  the  3-phase  into 
a-phase.  A  similar  claim  was  made  for  w  values  in  the  a-phase  near  w1 
and  w  values  in  the  3-phase  near  w^  as  to  the  existence  of  a  slowly 
moving  surface  propagating  from  the  a-phase  into  the  3-phase.  In  this 
section  we  will  give  an  argument  based  on  local  asymptotics  which  provides 
some  justification  for  those  claims. 

Let  us  consider  the  first  case  mentioned  above.  When  in  the  case 

when  <  Y  and  slowly  increasing.  We  assume  there  is  a  singular  surface 
located  at  point  x  =>  yQ  >  0  at  time  t  -  tg  where  w_,  u_  (w  <  w^)  and 
w+,  u+(w+  >  w2*  are  the  of  w,  u  as  x  yQ  from  below  and  above 

respectively.  As  we  are  interested  in  the  resolution  of  the  discontinuity 
in  w  in  am  immediate  neighborhood  of  x  -  yQ,  t  *■  tQ  we  introduce  stretched 
coordinates 


where  e  >  0  and  small. 

If  we  set  U(X,T;E)  *  u(x,t),  W(X,T>e)  ■  w(x,t)  we  find  via  the  chain 
rule  that  (1.8)  implies 


We  look  for  an  asymptotic  expansion  of  U,  W  in  £  of  the  form 


U(X,T»C)  -  Uq(X,T)  +  EU t ( X , T )  +  •• 
W(X,T?e)  -  Wq(X,T)  +  EW 1 ( X , T )  +  •• 
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Formal  matching  orders  of  e  in  (4.1)  shows 


“<V'X 


(4.2) 


If  U(X,T;e)  and  W(X,T;£)  are  to  approximate  the  behavior  of  u,  w  in 
a  small  neighborhood  of  (tg/Yg)  we  shall  need 


Ug(X,0)  -  u_,  WQ(X,0)  -  w_,  X  <  0 
UQ(X,0)  -  u+,  WQ(X,0)  «  w+,  X  >  0 


(4.3) 


Thus  (4.2)/  (4.3)  represents  a  Riemann  initial  value  problem  for  Ug,  Wg. 

A  solution  (4.2),  (4.3)  possessing  admissible  singular  surfaces  may  be 
constructed  in  a  manner  similar  to  a  method  given  by  Lax  [18]  for  hyperbolic 
conservation  laws.  Basically  the  idea  is  to  connect  u_,  w_  and  u+,  w+  to 
nearby  values  by  shock  or  rarefaction  waves.  The  nearby  states  can  be  then 
joined  by  a  slowly  moving  admissible  inter-phase  singular  surface  if  w+  is 
near  w2,  w+  >  w2  and  w_  is  near  w^,  w_  <  wq.  A  similar  construction  can 
be  done  for  w_  near  w^ ,  w_  <  w^ ,  w+  near  w^,  w+  >  w^.  These 
constructions  may  be  found  in  the  paper  of  James  [19].  Thus  for  >  0, 

sufficiently  small,  and 

-e,  <  »_  -  «,  <  0.  0  <  -  *2  <  I«+  -  “.I  <  *, 

(-c,  <».-«,<  0,  0  <  *+  -  »6  <  I»t  *  “-1  <  *l> 

(4.2),  (4.3)  possesses  a  solution  possessing  admissible  shocks  and  an  inter¬ 
phase  singular  surface  moving  with  speed  T  <  0  (T  >  0). 

We  note  similar  constructions  have  been  given  in  the  forthcoming  paper  of 
M.  Shearer  [20]. 

The  solution  Ug,  WQ  provides  an  asymptotic  estimate  of  the  solution  of 
(1.8)  in  the  neighborhood  of  the  singular  surface. 


5.  Discussion 


We  have  given  a  theory  of  visco-elastic  fluid  flow  which  exhibits 
hysteresis.  Our  theory  also  shows  the  fluid  exhibits  a  type  of  shape 
memory.  The  theory  is  qualitatively  similar  with  respect  to  the  hysterisis 
loop  as  the  experiments  described  in  [ 1 ] .  We  do  not  know  whether  the  fluids 
described  in  [1]  will  also  exhibit  shape  memory. 
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